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Abstract: It is generally held that neutrinos with superluminal velocity will lose their 
energy spontaneously by radiating electron-positron pairs, similar to bremsstrahlung 
process. Recently, this process was closely studied for neutrinos whose energy is roughly 
proportional to their momentum. Confronted with an increasing amount of superlumi- 
nal neutrino models, it is urgent to calculate the same process for general dispersion 
relations. The calculation is performed in this paper, without resorting to any nontrivial 
frame such as the effective "rest frame" . 
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1 Introduction 

The OPERA experiment stirred the physics community recently with its astonishing 
result that neutrinos in this experiment travel apparently faster than light at a high 
confidence level [1]. If not attributed to systematic errors in the measurement, this 
result would imply the violation of special relativity. Subsequently, it has inspired a 
lot of speculation. 1 

However, as quickly claimed in [21, 22], several high-energy processes disfavor 
the superluminal interpretation of the OPERA data. An outstanding example is the 
bremsstrahlung-like process 

^-^^ + e + + e", (1.1) 

where electron-positron pairs are radiated and hence neutrinos lose their energy effi- 
ciently To explicitly demonstrate this point, the authors of [21] assumed a special 
dispersion relation roughly of the form E = v u p, where v v is a constant greater than 
light velocity. A similar assumption was also taken in [22]. 

1 As a partial list, see [2-20] and references therein for speculation on this issue from various aspects. 
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On the other hand, a dispersion relation of the form E = v u p is too oversimplified 
to accommodate more observational data of neutrino velocity, as summarized and an- 
alyzed in [23]. Therefore, it is urgent to extend the calculation of [21, 22] to general 
dispersion relations. The present paper is devoted to such a calculation. Our results 
would help to rule out more phenomenological models and hunt for viable models. 

In the absence of Lorentz invariance, the calculation is complicated even for the 
special dispersion relation E = v u p, so an effective "mass" was assigned to neutrinos 
and a "rest frame" was employed in [21]. Our calculation does not resort to such a 
nontrivial reference frame. Or more explicitly, one may interpret our work as a direct 
calculation in the laboratory frame. Applied to the above dispersion relation, our result 
provides a crosscheck for the results in [21, 22]. 

The paper is organized as follows. Section 2 collects the basic assumptions and 
conventions of notation in this work. Kinematically there is a threshold energy for 
the process (1.1). We discuss the dependence of this threshold on dispersion relations 
in section 3. Section 4 is the main part of our paper, where we calculate the "decay 
width" of superluminal neutrino via (1.1). Details and techniques are presented clearly. 
To check and apply our general results in sections 3 and 4, we work out some specific 
examples in section 5 with given dispersion relations. The results are consistent with 
[21] quantitatively and [26] qualitatively. Section 6 concludes this paper. 

2 Basic assumptions and notation conventions 

Throughout the paper, we assume 

1. The ordinary conservation law for energy and momentum is intact. In other 
words, the time and space translations are exact symmetries in the working frame. 
A case study for violating this assumption can be found in [24]. 

2. The space is Euclidean and isotropic. Thereby we can work in a spherical coor- 
dinate system and define the magnitude of momentum as p = \p\. 

3. In the relevant energy range, the dispersion relation of electron and positron is 
well characterized by E 2 = p 2 + vr? e . This assumption is in accordance with 
experiments to date. 

4. The dispersion relation of neutrino is either E = E(p) or p = p(E). Here E(p) 
and p(E) are arbitrary devisable functions. They could be non-monotonic and 
may involve some parameters such as mass, etc. 
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In section 3, when deriving the threshold energy of process (1.1), we make one 
more assumption: 

• The neutrino's dispersion relation reduces to E 2 = p 2 + ml at very low energies, 
typically lower than the threshold energy by a factor O(10~ 4 ). See details in 
section 3. 

In section 4, this assumption is replaced by two other assumptions: 

• The squared amplitude (4.2) is the same as that in standard model. We will 
comment on possible loopholes of this assumption in section 6. But to alleviate 
the complexity of our calculation, we should make such an assumption at the 
moment. 

• The masses of electron and positron are neglected. This assumption seems to be 
reasonable, because when we study the decay width of high-energy neutrino which 
is practically around or above GeV scale, the phase space should be dominated 
by high-energy particles in principle. 

Let us clarify the conventions of notation by writing (1.1) in the form 

v(p) -> v(p') + e+{k') + e~{k), (2.1) 

where we have specified the notation of momentum for each particle. One must be 
careful of the notations of momenta. Taking p for instance, sometimes it stands for 
the four-vector, but sometimes it stands for the magnitude of three-vector p. In our 
equations, the four-vectors appear usually together with a dot, indicating the inner 
product with Lorentz signature diag(+, — , — , — ). For example, p ■ p = E 2 — \p\ 2 but 
p 2 = \p\ 2 - So there are no confusions if one is careful. 

The velocity of light is a constant, so we set it to 1 throughout this paper. 

3 Energy threshold 

The assumptions made in the previous section simplify the derivation of threshold 
energy for (1.1) considerably. 2 From the energy conservation relation 

E p = E p i + E k > + E k , (3.1) 

2 When our work was in preparation, ref. [25] appeared. The topics in this section and [25] overlap 
partly. For the completeness of our paper, we keep this section in its own form. Ref. [26], starting 
form different assumptions, also concerns a partly overlapped subject of this paper. It appeared more 
recently when we were polishing our work. 
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we observe that to minimize E p , one should lower Ep>, E k i and E k as much as possible. 
This can be achieved by deleting all of the transverse momentum components, leading 
to the reduced conservation law of momentum 

p = p ' + k' + k, (3.2) 

where p — \p\, p = \p'\ and so on are magnitude of momenta as our conventions. 

Substituting (3.1) and (3.2) into the dispersion relation of e + (k') and e~(k), we 
have 

(E p - E p ,f - 2E k (E p - E pl ) - (p - p'f + 2k(p - p') = 0. (3.3) 

In light of dispersion relations of v(p), v(p') and e~(k), this equation may be understood 
as an implicit function of E p , E p i and E k . Then we can extremize E p with respect to 
Epi and Ek, obtaining 

(E p -E k )-(p-k)^f = 0, 

(E p -Ep>)-(p-p')^ = 0, (3.4) 

where conditions dE p /dE p i = dE p jdE k = and dispersion relations dp'/dE p i = E p >/p', 
dk/dE k = E k /k are used. 

From eqs. (3.3) and (3.4), it is not hard to get 

2E k = Ep-E p i, 2k = p-p', ^p = Elf (3.5) 



p p' 



and subsequently 



E 2 - 2E p E p/ + ml = p 2 - 2pp' + 4m 2 e , E p , = J^-^ . (3.6) 



1 — 

v 



As a result, we find the threshold of (1.1) is given by 

(E 2 p -p 2 ) thr ={2m e + m v ) 2 . (3.7) 

This result is in accordance with [25]. To the leading order, it is also compatible with 
[21] where neutrino's dispersion relation is E p /p = v v . This will be shown in subsection 
5.2. As a concrete example, in subsection 5.3 we will utilize (3.7) to get the decay 
threshold for a toy model of mass-dependent Lorentz violation [23]. 

Note that when writing down (3.4) and (3.6), we have assumed the dispersion 
relation E 2 , = p' 2 + m 2 u of v(p') at low energies. But the dispersion relation of 
v(p) is irrelevant throughout the derivation. From eqs. (3.6) and (3.7) we can see 



-4- 



[E v i I -Ep)thr. = m v /(2m e + m u ) ~ C(10 -4 ). That means we have assumed the disper- 
sion relation E 2 , = p' 2 + m 2 for neutrinos at an energy lower than the threshold energy 
by (9(10 -4 ). This assumption is natural and consistent with the observational data of 
neutrino velocity summarized in [23]. Replacing this assumption with other dispersion 
relations, one may also restart from (3.3) and follow our method to derive the threshold 
energy. 3 



4 Decay width 



We are interested in the following process: u(p) — > u{p' 
the neutral current of this process, we have 4 



^MM* = 128G 2 F [(p ■ k')(k ■ p')(~ 



sin 



e + (k') + e (k). Considering 
2m2 1 (p-k)(k' -p') sin 4 6)}. (4.1) 



spin 



To calculate the decay width, we will integrate over p', k and k', hence we can make 
use of the symmetry between k and k! to write the squared amplitude as 



MM* = 128G 2 F (p ■ k'){k ■ p') 



spin 



- + sin 2 6 W 



+ sin 4 9 W 



(4.2) 



The decay width is formally given by 
1 f d 3 p' f d 3 k 



2E P J (2vr) 3 2£p/ J (2n) 3 2E, 



d 3 k> 



8G 2 F 



(2tt) 5 E p \A 



— sin 9 W + 2 sin 4 9 



(2n) 3 2E k/ 2 

d 3 p' f d 3 k 
Ek 



-\M\ 2 {2n) A 5 4 {p — p' — k' 
d 3 k' 



k) 



(4.3) 



E p i 



2Ei 



(p ■ k')(k ■ p')5 4 (p - p' -k' -k). 



Here 8w is the Weinberg angle. 

In the main part of this section, we will focus on the calculation of integral by 
temporarily forgetting the overall coefficient. Because we are interested in high-energy 
neutrino decay, the masses of electron and positron will be neglected in our calculation. 



3 This method is valid if E p has local minima as an implicit function (3.3) of £y and Ek- The 
situation will be more complicated if the configuration of (3.3) does not have a local minimum. 

4 In the first version of our manuscript, the squared amplitude is incomplete. Here we corrected this 
error and included the last term in (4.2). The difference does not affect most of our calculations. It 
only modifies an overall factor in the decay width. We are grateful to Zhaohuan Yu, Fedor Bezrukov 
and Evslin Jarah for communications on this point. 
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Then the integral in (4.3) can be briefly rewritten as 
d 3 p' f d 3 k f d 3 k' 



abb. 



? -( P .k')(k.p')5\p-p'-k'-k) 
f ^ f d * k > 6 {k > . k >) \ klo>o{p . k ' ){k . p > )5 * (p _ p > _ k > _ k) 



E p > J E k 
d 3 p' f d 3 k 



Ep> J E k 



(4.4) 



[p-(p-p' - k)] (k ■ p') 5 ((p -p -k) ■ (p-p - k)) \ Ep „ Ep ,-k>o- 



Making use of relations k ■ k = k' ■ k' = and 
k ■ p' = k ■ (p - k - k') = k ■ p - {k + k>) ' {k + k>) 
we obtain 

d 3 p' f d 3 k 



k ■ p 



(p - p') ■ (p - p') 



(4.5) 



abb. 



E p i 



[p ■ (p — p) — p ■ k] 



k ■ p 



(P - P') ■ (P - P') 



x5 ((p — p') ■ (p — p') — 2k-(p — p')) 
d 3 p 



Epi 



x 



E p E k -p-k 



E p -E p ,-k>0 

kdk sin 6\d6\d^)\ p ■ (p — p') — E p E k + p ■ k 
(p - p') ■ (p - p') 



xS((p- p') ■ (p - p') - 2E k (E p - E p ,) + 2k ■ {p-p') 



(4.6) 



E p -E p ,-k>0 



Here 9\ is defined as the angle between k and p — p' . We will define 82 as the angle 
between p and p'. The relative directions and angles between the relevant momentum 
vectors are depicted in figure 1. It is convenient to express the inner product of p and 
k in terms of the new coordinates, 



p • k = pk 



cos 6*i (p — p' cos 62) + p' sin 9\ cos (pi sin 82 
\p — p'\ 



(4.7) 
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Figure 1. (color online). The relative directions of vectors p, p', p — p' and k in three spatial 
dimensions. The zenith angle 9\ is defined as the angle between k and p — p', and #2 as the 
angle between p between p' . As shown in the picture, axes e 2 and e z > coincide with directions 
of p and p — p 1 respectively. Projecting k on the plane perpendicular to e z r, we can define 
one azimuth angle ip±. Similarly, the other azimuth angle ip2 can be defined by projection of 
p' on the plane perpendicular to e z . 



In terms of the new coordinates, the integration takes the form 
d?p' 



abb. 



kdk sin d\dd\d<j)\ 



P- (P-P') ~ E P E k 



-pk 



cos 6\ (p — p' cos 62) + p' sin Q\ cos tpi sin 62 



E p E k - pk 



\p — p'\ 

cos 6*i (p — p' cos 6 2 ) + p 1 sin 61 cos ipi sin 9 2 (p — p') ■ (p — p') 



\p — p'\ 



xS{{p- p') ■ (p - p') - 2E k (E p - E pl ) + 2k\p- p'\ cos^i 
d 3 p 



E p -E v ,-k>Q 



E v i 



kdkdifi i dx 



P-{p-p) - E p E k 



+pk 



x(p — p' cos 62) + p'yl — x 2 cosifi sin 9 2 



\p — p'\ 



E p E k - pk 



x(p — p' cos^) + p'V^ — x 2 cos <f 1 sin 9 2 (p — p') ■ (p — p') 



\p — p'\ 



x5 ((p - p) ■ (p - p) - 1E k (E v - E pl ) + 2k\p- p'\x 



(4.8) 



E p -E p ,~k>0 
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With the newly introduced variable x = cos#i G [—1, 1], we note that the function 
f(x) = {p- p') ■ (p - p') - 2k{E p - Epf) + 2k\p- p'\x (4.9) 

has a single root 

x = 2k (E P - E p ,) - (p -p>) ■ (p - p>) 
2k\p — p'\ 

and its first-order derivative 

f\x) = 2k\p-p'\. (4.11) 
Therefore we can integrate x out of the delta function and quickly obtain 

d 3 p 



abb. 



E p i 



kdkdipi 



p-(p-p') - kE p 



+kp 



x {jp — p' cos 9 2 ) +p' a/1 — x\ cos cpi sin 9 2 



\p — p'\ 



x [E p k — pk 



xo(p — p' cos 9 2 ) + p' y l — Xq cos if i sin 9 2 (p — p') ■ (p — p') 



\p — p'\ 



x- 



Ep-E pl -k>0 



dki 2vr 



p • (p — p') — kE p + p 



-P- 



2k\p — p'\ 

d 3 p' 
2E p ,\p-p'\ 

kx (p — p' cos 9 2 ) {p — p') • (p — p') 
\p — P'\ 2 



71 



kx (p — p' cos 9 2 ) 
\p — p'\ 
(pp') 2 (k 2 - k 2 x 2 ) sin 2 9 2 
\p — p'\ 2 



E p k 



(4.12) 

E p -E pf -k>0 



One may check that the integrand of (4.12) is independent of <p 2 . Its dependence on 
k is quite simple. Its domain of integration is determined by — 1 < x < 1, E p —E p i — k > 
0, or namely 



(p - p') ■ (p - p') 
2{Ep-E pl ) + 2\p-p'\ 



< k < 



(p - p') ■ (p - pQ 

2(Ep-E p ,)-2\p-p>\ 



E p - Ep, > \p-p'\. 



(4.13) 



So we can integrate over variables k and tp 2 straightforwardly. The calculation is a 
little tedious, yielding 

p' 2 dp'dy 



8G 2 F 



sin 9 W + 2 sin 9 W ] n 




6E P , 



3E p Ep, 



{2-nfEp \A 

+ (-QE 2 p + 2p 2 - 3pp'y)E 2 + (3E 3 p - 3E p p 2 - 3E p p' 2 + 8E p pp'y)E p 



+3pyp' i + (2E 2 p - 2p z - ApY)p' 2 + (3p A y - 3pE I p y)p 



(4.14) 
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The domain of integration is 

p 2 + p' 2 - Ipp'y < (E p - E p ,f, (4.15) 
-1 < y < 1. (4.16) 

In the above, we employed notation y = cos 9 2 and turned on the overall factor in front 
of the integral as given by (4.3). 

The expression (4.14) for decay width is one of our main results. In section 5, we 
will apply it to several examples with explicit dispersion relations. Before proceed, let 
us make some comments on (4.14). 

First, we would like to show that the decay width (4.14) is positive-definite, as one 
should have anticipated from its definition (4.3). For this purpose, it is enough to focus 
on the integrand inside the square brackets, which can be transformed to 

[■■■]= 3{E p E p , - pp'y) [(E p - E pl f - (p 2 + p' 2 - 2pp'y)} + 2(E p p' - E p ,p) 2 

+2pp'(l - y) [2E P E P , - pp'{\ + y)\ . (4.17) 

Taking account of inequalities (4.15) and (4.16), this expression is nonnegative if 
E p E p i > pp' . This condition is well-satisfied if E p — p > 0, E p > — p' > at lower 
energy and dE p /dp > 1, dE p i/dp' > 1 in the energy region of superluminal neutrino. 

Second, we note that (4.15) puts a lower limit of integration y > [p 2 + p' 2 — (E p — 
Ep/) 2 } / (2pp') . In most situations, we have E p — E p i < p + p', then this limit is more 
stringent than y > —1, and hence the domain of integration is simply [p 2 +p' 2 — (E p — 
E p i) 2 ]/ {2pp') < y < 1, leading to the reduced decay width 

8Gf (I . 2 . 4 \ 2 [p' 2 dp> 1 r 6 

r Ep - Ep ,< P+ ^ = [~ A ~ sm 9 W + 2sm 9 W J n J — — {5E p 

-3 (15E 2 + 5p 2 - 3p' 2 ) El + 8E P , (lOE 2 - Qp' 2 + 9pp') E* p 

-3E 2 [15E% - 6 (3p 2 - 8p'p + 3p' 2 ) E 2 p , - (p - p' f (5p 2 + Wp'p - 3p' 2 )] 

-24E p E p ,p [{2p - 3p')E 2 + 3(p - p'fp'] + hE % p , + 3E A pl (3p 2 - hp' 2 ) 

-3E 2 pl (p - p') 2 (3p 2 - lOp'p - hp' 2 ) - 5(p - p'f (p 2 + Ap'p + p' 2 ) }.(4.18) 
5 Examples 

In sections above, we have derived the kinematical threshold and "decay width" of 
the bremsstrahlung-like process (1.1) for general dispersion relations of neutrino. This 
was done under the assumptions made in section 2. To check our main results (3.7) 
and (4.14), we will apply them to muon decay process in subsection 5.1 and to Cohen- 
Glashow model in subsection 5.2. As further applications, we will use them to study 
some other models in subsections 5.3, 5.4, 5.5. 
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5.1 Muon decay 

Process (1.1) can be regarded as a three-body decay process by weak interaction. It 
is analogous to muon decay in standard model. So the basic test of our result (4.14) 
is comparison with /i(p) — > v^{p') + v e (k') + e{k) by neglecting particle masses in the 
final states and replacing neutrino with muon in the initial state. Setting E p i = p', 
E 2 = p 2 + m 2 , we work out (4.14) directly 



G 2 m 6 ( 1 

sin 2 # VK + 2sin 4 fl w ) . (5.1) 



1927r 3 .Ep \4 

It is different from the muon decay width by a factor (1/4 — sin 2 9yy + 2 sin 4 9\y)- This 
factor should be replaced by 1 if we incorporate the charged current. So our result 
exactly passes the muon decay test. 

5.2 Cohen-Glashow model 

The second test is to recover the result in [21]. For this purpose, we set E 2 = p 2 (l +8), 
E 2 , =p' 2 (l + 5) in (4.14) and get 

4 G 2 E 5 5 3 f 1 \ 
F = - x (_ _ sin2 9w + 2 sin 4 9 W J . (5.2) 

The decay width of [21] can be numerically reproduced 5 by taking sin 2 6w — 1/4. 

Another main result of this paper is the threshold (3.7). Applying this threshold 
condition to the model of [21], we find 

2m P + m v 2m P 
E thl . = * ^ (5.3) 

the same as the threshold in [21]. 

5.3 Mass-dependent Lorentz violation 

In ref. [23], two of the authors proposed the mass-dependent Lorentz violation scenario 
to explain the observed neutrino velocity as a function of energy. In this scenario, the 
mass-energy relation of neutrino has the form 

l-v 2 = A -/(A), A = m 2 /£ 2 (5.4) 



5 This corrects the wrong claim in the first version of our manuscript, because the modification of 
squared amplitude changes the overall factor in the decay width. 
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2m e / GeV 



2m e /GeV 



Vl-p 2 /£ 2 y/l-p 2 /E 2 




E/GeV 



0.001 0.1 10 1000 0.30 0.35 0.40 0.45 0.50 0.55 0.60 



Figure 2. (color online). Numerical solution of the threshold for a toy model of mass- 
dependent Lorentz violation. The solid blue line depicts function 2m e (l — p 2 / E 2 )~ 1 / 2 . It 
crosses the dashed black line at the threshold energy, as highlighted by a black dot in the right 
graph. The electron/positron mass is set to 0.5 MeV. The toy model and other parameter 
values are chosen to be the same as in [23] . 



where /(A) is a mo del- dependent function. The new function /(A) is useful phenomeno- 
logically, because we can get its information directly from experiments which con- 
strain neutrino velocity as a function of energy. With the definition of group velocity 
v = dE/dp, relation (5.4) can be taken as a differential equation and integrated into 
dispersion relation 



i m 2 iff m? \ 

E2 "I - J y& J 

A concrete toy model of mass-dependent Lorentz violation was devised in [23], well 
fitting observational data of neutrino velocity. For the toy model and parameters given 
in [23], we combine this relation with eq. (3.7), and numerically get the threshold 
energy at about 0.5 GeV. This is illustrated in figure 2. This threshold is higher than 
that of the Cohen-Glashow model [21]. 

Because the dispersion relation takes a very complicated form, it is difficult to work 
out the decay width in this toy model, even numerically. As an alternative, we will 
deal with a simplified model in subsection 5.4. 

5.4 Velocity of step form 

In the toy model of ref. [23], the dependence of velocity on energy looks like a delta 
function, and well explains the observational data of neutrino velocity. But it is very 
difficult to work out the decay width (4.14) for that model. As an alternative, let us 
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r/r CG 

1.0 r 




E/GeV 



Figure 3. (color online). The dependence of decay width (4.14) on neutrino energy Ep in 
model (5.6). We have fixed 5 = 2 x 1CT 5 , E c = 5 GeV for the blue points and E c = 17 GeV 
for the red points. 



study a model in which the neutrino's velocity depends on energy in a step form 

dp/dE = l+( 1 - 1 ) H{E-E C ). (5.6) 



VT 

Here H(x) is a Heaviside unit step function and E c is an critical energy. 

Fixing 5 = 2 x 10~ 5 , we numerically computed the decay width and got the results 
in figure 3. From the figure we can see the value of decay width becomes closer and 
closer to (5.2) as E increases. The figure also tells us that the decay width gets larger 
if neutrino is faster than light in a wider energy range. Reversing the logic, the decay 
width can be suppressed by narrowing the energy range in which the neutrino is faster 
than light. Perhaps this is realizable torturously if the dependence of neutrino velocity 
on energy takes a comb-like form. 

5.5 Horava-Lifshitz model 

Motivated by Horava-Lifshitz theories, ref. [26] has studied the dispersion relation for 
neutrinos of the form 

E 2 =p 2 + m 2 + r ] 'p 2 + ^-. (5.7) 

Here the mass of neutrino m is negligible. When the r( correction dominates, this 
model reduces to the Cohen-Glashow model. 

When the 77 correction dominates, dispersion relation (5.7) becomes 

E2 =p 2 + w ( 5 - 8 > 

with energy-dependent velocity. For this form of dispersion relation, we can calculate 
the decay width with (4.14). To the leading order of 77, it is 
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1665 fE 2 \ 3 G 2 F Elrf (1 
r = — x U x 1 P ~- sin 2 0^ + 2 sin 4 fry . (5.9) 



2002 \M 2 J 192vr 3 V 4 

Remembering that the neutrino velocity is v 2 — 1 ~ 3rjp 2 /M 2 , it is convenient to take 
the notation 5 = 3r\E 2 jM 2 and rewrite decay width (5.10) as 

„ 185 G 2 E 5 J S (\ 



x 



2 E b d 6 ( 1 \ 

/ _ sin 2 ^ + 2 gin 4 Q\ (510) 



6006 192tt 3 

At energies relevant to the OPERA experiment, such a decay width is one or two 
orders smaller than that of the Cohen-Glashow model. This ratio of suppression is 
consistent with the results of ref. [26]. A naive application of (3.7) to (5.8) yields 
?7P^ hr /M 2 = (2m e + m u ) 2 , which gives the threshold energy in leading order 

£th, * -yi ■ (5.11) 
It is unsafe to take -E t hr. — 2m e \/3/5, because here S is energy- dependent. 



6 Conclusion 

In this paper, under the assumptions enumerated in section 2, we studied the kinematic 
threshold and decay width of superluminal neutrinos for the bremsstrahlung-like pro- 
cess (1.1). This was done for general dispersion relations of neutrino, without resorting 
to any nontrivial frame such as the effective "rest frame". The main results are rep- 
resented by eqs. (3.7) and (4.14). Our results confirmed and generalized the previous 
results in [21, 22]. 

Before concluding this paper, we would like to make some relevant remarks on the 
assumption of squared amplitude, which leaves a loose end for future investigation. As 
has been emphasized in section 2, when calculating the decay width, we assumed that 
the squared amplitude is the same as that in standard model. Strictly speaking, this is 
not always a consistent assumption when general dispersion relations are involved. In 
general, dispersion relations will enter into both kinematics and dynamics of particle 
physics. However, without an assumption on the squared amplitude, we cannot do 
any calculation about decay width. At the same time, since the deviation of neutrino 
dispersion relation is not too far from special relativity, we expect that the deformed 
amplitude should not deviate significantly from the standard model. Therefore, our 
assumption (4.2) is not only necessary but also natural to some extent. We thus expect 
our result provides a good estimation of decay width in order of magnitude. Of course, 
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further investigation is required to confirm this expectation and improve the present 
situation. See ref. [27] for a recent progress along this direction. 

Another interesting project is employing our general results to rule out more phe- 
nomenological models and hunt for viable models, as shown by some examples in section 
5. We feel this project will be challenging but rewarding, given the importance of special 
relativity in modern physics. 
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